428 



The Sf ability of Fluid Motion. 
By T. H. Havelock, F.RS. 

(Eeceived January 31, 1921.) 

1. The following notes on the stability of fluid motion arose from a desire 
to use the energy method, introduced by Reynolds and modified by Orr, as a 
measure of the comparative degree of stability of various types of flow under 
different boundary conditions. A few examples are worked out to illustrate 
this point of view : in § 5 a case which resembles the flow of a stream with a 
free surface ; in § 7 flow which approximates to a uniform stream between 
fixed walls without slipping at the walls; in §§6,8 motion with other 
boundary conditions. Before proceeding to these, it seems desirable to give. a 
short account of the method in the form in which it is used later, together 
with some remarks on its relation to the classical method of small vibrations. 

2. We shall consider only two-dimensional motion of an incompressible 
viscous fluid limited by the planes y = ±a. Let the steady state under an 
assigned forcive and given boundary conditions be specified by a velocity, U, 
parallel to the axis of x. Let the disturbed state have velocity components 
(TJ + u, ^)and let the additional pressure be p. Then, by taking the difference 
of the two sets of hydrodynamical equations for the two states and neglecting 
squares and products of the additional velocities, we have 

du , TT du , 3U 1 dp , a 9 
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together with the equation of continuity. 

It is convenient to introduce non-dimensional variables given by 

x = a% ; y ~ arj : ar = XJt ; 

where U is the mean velocity over the cross-section in the steady state. 
Further, we write TJU instead of U, and take the current function of the 
additional velocity to be XJa^jr. Eliminating^ from the two equations (1), we 
obtain 

E (^ + |)^- EU "t ==2vV ' (2) 

where U /; is written for d 2 V/d?) 2 , and R is Reynolds' number 2aU/p. There 
are in addition the appropriate boundary conditions for the disturbing 
function, ty, The classical method of examining tfie stability of a given 
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distribution U consists in assuming a solution of (2) of the form 
exp. {i(wr +pl~)} fn(v)- ^ or an y arbitrary real value of p, the corresponding 
possible forms of f n (v) an d values of n are found from (2) together with the 
boundary conditions. The distribution TI may be said to be thoroughly 
stable if every possible value of n has a positive imaginary part, and if this 
holds for all positive values of E. 

The usual boundary conditions, which we shall assume in the first place, 
are u = 0, v = 0, or 

^ = 0; df/d v = 0: v = ±1. (3) 

From the work of Kelvin, Kayleigh, Orr, Hopf, and others, it may be taken 
that the simple shearing motion, U = 1-f ?/, is thoroughly stable in this sense ; 
and probably a similar conclusion holds for motion under a constant force or 
pressure gradient, namely U = f (1 — if). 

There are various possible explanations of the well-known divergence 
between these results and the behaviour of actual fluids. In the first place, it 
is obvious that the physical properties, whether of the fluid or of the walls, are 
inadequately specified in the mathematical statement of the problem. But, 
apart from this, the disturbances have been supposed small, and second order 
terms neglected. Again, in a system of this type, a disturbance may be 
small initially and may converge ultimately to zero, but may be very large at 
intermediate times, and may thus give rise to practical instability. 

The energy method of Reynolds is in a different category from these in 
that it takes the mathematical problem as it stands and does not necessarily 
involve the actual magnitude of the disturbance ; in fact, it forms a new 
criterion or measure of degree of stability. The energy of the disturbance 
being defined by _ rfr /^\2 fiW) 

we have from (2) and (3), after integrating by parts, 



6 i5 _ jL„r?2 
dt ~ 2fX 
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(5) 



Here dE/dt means the rate of increase of E in a region whose end boundaries 
move with the steady velocity U. We may replace this by dE/dt for a region 
with fixed ends, and we shall then have additional terms on the right of (5) 
denoting flux of energy across these ends. The latter terms may be omitted 
under conditions which cover the usual cases : namely, either the disturbance 
is periodic in f , or it is limited or localised so that ^ and its derivatives 
converge sufficiently rapidly to zero for £ = ± oo . "We shall assume such 
conditions to hold in what follows, and references to boundary conditions 
mean those which hold at the planes tj = +1. 
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Eeynolds' method of using (5) to determine a criterion of stability 
consisted in assuming a suitable form for i/r and finding the least value of K 
for which the right-hand side of (5) is zero. It is usually stated that this 
method assumes turbulent motion to be already in existence, and it then 
gives a criterion to show whether the turbulence is increasing or decreasing 
momentarily; but this is somewhat misleading without defining what is 
meant by turbulent motion. Equation (5), as stated above, applies to any 
small arbitrary disturbance, neglecting terms of the second order, as in the 
ordinary method of small vibrations; further, U is a laminar fluid motion 
satisfying the usual hydrodynamical equations under the given conditions. 

On the other hand, Eeynolds defined XJ as the mean velocity at each 
point, taken over a small region or during a short time, and this principal or 
mean motion need not satisfy the ordinary equations. The extra velocities 
u and v then play a double part, in that they specify the disturbance, and at 
the same time give a measure of the turbulence ; they must satisfy certain 
conditions as to their mean values, and then equation (5) holds in the same 
form when mean values are used. However, in applying it to find the 
criterion for flow under a constant pressure gradient, Eeynolds, and Sharpe 
following him, did, in fact, take U to be the usual form, C (a 2 —y 2 ), for steady 
laminar flow. But in turbulent flow, although the variations of velocity at 
any point are small, yet they may cause the gradient of the mean velocity to 
differ appreciably from its value in laminar flow, as is obvious from a com- 
parison of the curves of distribution of velocity across a pipe in regular and 
in turbulent flow. 

However, it is unnecessary to dwell on this distinction, as it has been 
pointed out clearly by Lorentz* and other writers ; further, we shall 
consider here only small disturbances. 

3. Under these circumstances, the energy method has been given a precise 
and definite meaning by Orrf from the following considerations :— 

If the right-hand side of (5) is positive, the energy of the disturbance is 
momentarily increasing. But, for a given velocity distribution, 17, it may be 
impossible to find any function, yjr, satisfying the boundary conditions, such 
that that expression is positive, unless E exceeds a certain value. If such be 
the case, this least value of E is a critical value of definite significance. The 
corresponding critical disturbance is found by taking the variation of the 

eqimtl0n E f JU'II -^ dgch,-2 \[(v 2 f) 2 d£d v = 0, (6) 

subject to SE = 0. 

* H. A. Lorentz, ' Abhandlungen fiber Theor. Phys., 3 vol. 1, p. 43. 
t W. Mc.F. Orr, < Proc. Roy. Irish Acad., 1 vol. 27, p. 9 (1907). 
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Carrying out the variation, and using the boundary conditions (3), we 
obtain 

4 V V + 2KU , ^ + KU // ^ = 0. (7) 

C£ or) d£ 

To find the critical value of K, we assume first that £ occurs in ^asa 
factor exp. ipi*, and then solve (7) ; using the boundary conditions, we have 
an equation from which we can find the least value of R for a given value 
of p, and finally we take the minimum value of R with respect to p. 

The process has been expressed in a different form by Hamel.* Using the 
corresponding Green's function for the equation v 4 ^ = ®> tne equation (7) 
may be replaced by a linear integral equation for ty, of which the required 
value of E is the lowest characteristic number. 

Eeturning to equation (5), if dE/dt is positive for any assigned initial 
disturbance, it does not follow that the motion is unstable in the ordinary 
sense. But, if there exists an absolute minimum for R in the manner 
explained above, it follows that, when R is less than this value, dE/dt is 
negative for every initial disturbance, and must always remain negative. 
Thus the system has at least a much higher degree of stability for such 
values of R compared with those greater than the critical minimum. 
Obviously, this method does not produce any new information which is 
not implicit in the ordinary equations, such as equations (2) and (3) ; but it 
presents part of that information in a different form, so that the critical 
minimum of R may be used as a measure of the degree of stability of various 
distributions of velocity under different boundary conditions. 

4. It is convenient to classify the boundary conditions under which the 
energy equation (5) is valid. For this purpose we use an alternative form 
derived directly from equations (1), with the ordinary notation 

Pxx= - p+ 2„ Tx ; P* = /*(^t £ ); ft»=-* + 2/*£,. (8) 
We have 

dE Cr „ , „ , Cf dV 



dt 



{w(lpxx'+mpzy) + v(lpty + ^w)} ds—p I uv -r-dxdy 



+ 



du ?,v fdu dv\~] ■ 7 



where ds is a line element of the boundary and (I, m) the normal. 

We have specified .the conditions at the end boundaries, and we are con- 
cerned now with the planes y = ±a. It follows, that we get the energy 

* G. Hamel, 'Gott. Nachr,, Math. Phys. Klasse/ 1911, p. 261. 

2 h 2 
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equation (5), without any surface integrals expressing transfer of energy across 
the boundaries, with the following combinations : 

(i)u=0, v=zO; (ii) u = 0,^ = 0; (iii) v—0 t p^ = 0; (iv) p xy = 0, p yy = 0. 

We may also verify that, under these conditions, the variation of (6) leads 
to the same differential equation (7). 

5. Most of the fluid motions whose stability has been examined, come under 
case (i) of the above. A different case of special interest is a stream with a 
free upper surface, the conditions at the upper surface being as in (iv). These 
conditions, however, do not lead to simple expressions in terms of the 
disturbing function, i|r ; moreover it is not permissible to regard the upper 
free surface as rigorously plane. We therefore, following Kelvin,* replace 
the problem by one which is very nearly the same but is more easily 
specified; it may be described as a broad river flowing over a perfectly 
smooth inclined plane bed, the upper surface being fitted by a parallel plane 
cover moving with the water in contact with it. The conditions at the 
upper surface then come under case (iii) of the previous section. 

We take the origin in the ujDper surface in this case, so that a is the depth 
of the stream and E is aJJ/v. The steady state is given by 

U = #(l-V). (10) 

Using this in (7) and assuming ty to be proportional to <f^, the differential 
equation becomes 

£-lJ f -,(2J± + +)=0, (11) 

where a = prj, and k = 3iH/2p 8 . 

The boundary conditions are u = 0, v = at the bed of the stream, and 
v = 0, p xy = at the upper surface ; these reduce to 

ty = 0, d 2 yfr/dcn 2 = ; a = ; 

i/r = 0, dtyjdoL = ; a = p. (12) 

Equation (11) was solved by Orr for flow between two fixed planes with u 
and v zero at both boundaries, and it was found necessary to consider only 
solutions in even powers of a. We shall require here the corresponding- 
solutions in odd powers. Writing a solution in the form 

^r = 2 A n u n /n I 
we have the sequence relation 

A M+4 — 2A n+2 -f{l— (2n+l)k}An = 0. , (13) 

* Kelvin, ' Math, and Phys. Papers,' vol. 4, p. 330. 
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Denoting by ^Jr , ^i, yfr 2 , ^3 the solutions beginning with 1, ot, a 2 y a s 
respectively, it follows from (12) that the boundary conditions lead to 

TJndy]rsfdu--'tysdy}rifdoi=z 0. (14) 

where a has to be replaced by p. 

Calculating the coefficients far enough to give sufficient accuracy for our 
purpose, we have 

fi = <% + 2aV3J + (3 + 3ft)* 5 /5i+(4+20ft)*y7! 

+ (5 + 70ft + 33ft 2 )*y9! + (6 + 180ft + 366ft 2 )a u /lli 

+ (7 + 385ft + 2029F + 627F)a 13 /13l-f(8 + 728A + 7832F + 9672F)a 15 /15! 

+ (9 + 1260ft + 24030F + 73500ft 3 + 16929ft 4 )* 1 7l7 ! + ..*, 
^3= ^/3! + 2a 5 /5l+(3 + 7^)a77i+(4+36^)a 9 /9! 

+ (5 + 110ft + 105ft 2 ) a n /ll ! + (6 + 260ft + 894ft 2 ) « ia /13! 
+ (7 + 525ft + 4213ft 2 + 2415ft 3 ) a 15 /15 ! + (8 + 952*+ 14552ft 2 

+ 28968*»)««/17l+... 

Forming equation (14) we have 

2/3 I+8//5 I+32p*/7 ! + 128p«/9 !+(512 + 192F)//ll ! 
+(2048 + 2244F)y°/13 !+(8192 + 19456/i%> 12 /15 ! 
+(32768 + 139264P)/ 4 / 1 7 1 + (131072 + 901120P 

+ 129024^) p 16 /19 ! + ... = 0. (15) 

Only even powers of k appear in this equation, thus giving a check upon 
the arithmetic ; further, the terms independent of k may be summed. Taking 
the least root of (15) as an equation for & 2 , we have approximately 

•R2 _ . \ sinh2p--2p n fi . 

a «/192 , 2244» 2 , 19456p* , 139264»« , 901120» 8 , \* V ; 
^\11! 13! + 15! + 17! + 19! + **7 

Instead of forming an equation for the minimum value of K, it is simpler 
to find it by trial. We find, with sufficient accuracy, that it occurs near 
p 2 = 11, and then, approximately, 

E = 96. (17) 

The corresponding value, found by Orr, for flow under similar conditions 
but with a fixed plane at the upper surface, is 117. We conclude then that 
flow in an open canal has a lower degree of stability than flow between fixed 
planes. 

Turning to experimental results, the number usually quoted for flow 
through a tube is 2000 approximately. This was obtained chiefly from 
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experiments with smooth glass tubes ; a much lower number, of the order of 
400, has been found from metal tubes. The only available direct results for 
flow in an open stream appear to be those given *by Hopf,* who found E to 
be of the order of 300. These results agree in character with the theoretical 
calculations, which is all that could be expected. 

It is of interest to note that this appears to contradict a statement by 
Reynoldsf in one of his earlier papers. He classes separately circumstances 
conducive to steady motion and those conducive to unsteady motion : among 
the former a free surface, and in the latter solid bounding walls. However, 
this opinion seems to be based on visual observation of eddies caused by the 
wind beneath the oiled surface of water. "At a sufficient distance from 
the windward edge of an oil-calmed surface there are always eddies beneath 
the surface, even when the wind is light. . . . Without oil I was unable 
to perceive any indication of eddies." 

This introduces a different property of a boundary surface, namely, that of 
initiating disturbances. The mathematical statement ignores this property 
and specifies only control of the velocity functions : the disturbances are 
supposed to be initiated by some extraneous agency, and it is tacitly assumed 
that all types of disturbance are equally probable. It may be, for instance, 
that the theoretical results for flow through pipes should be compared with 
experiments on rough pipes rather than those with perfectly smooth walls. 
However, we may conclude that a solid boundary is conducive to stability in 
so far as it ensures that there is no slipping of the fluid in contact with it. 

6. In determining the minimum value of E from the differential equation (7), 
there are only two factors : the distribution of steady velocity, U, and the 
boundary conditions for the disturbance. The comparison in the previous 
section, between an open stream and flow between fixed walls, involved 
changes in both these factors. We may separate the effect of the boundary 
conditions by assuming the same value of U as in (10), but expressing the 
property of the supposed moving plane in contact with the upper surface by 
% = 0, v = 0, instead of by v = 0, f xy — 0. To anticipate the argument of 
the next sections, we should expect a value of E intermediate between 96 
and 117. 

We have the same equation (11) for ty, together with -xjr = 0, dyjr/du — 
at a = 0, and a = p. It follows that only the solutions i/r 2 and -\/r 3 are 
involved, and we have 

^jr 2 dtyzj 'da—^zd\l> r 2 / 'da = 0, (18) 

* L. Hopf, < Ann. der Phys.,' vol. 32, p. 777 (1910). 

+ O. Reynolds, l Scientific Papers/ vol. 2, pp. 157, 59. See also A. H. Gibson, 'Phil. 
Mag.,' vol. 25. p. 81 (1913). 
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when a = p. The series for yjr^ is given in § 5 ; also we have 

yfr 2 = a 2 /2!-f-2a 4 /4!4-(3 + 5^)a 6 /6! + (4 + 28/^)a 8 /8! 

+ (54-90^ + 65F)a 10 /10! + (6 + 220A+606P)a 12 /12! 

+ (7 + 455/^ + 3037F+1365F)a 1 Vl4!4-(8 + 840/^+10968P 

+ l7880F)a 16 /16 !+.... 
The boundary equation (18) leads to 
2/4! + 8p 2 /6! + 32^ 4 /8! + 128p 6 /10! + (512 + 280P)^ 8 /12! 
+ (2048 + 3136P) p 10 /14 ! + (8192 + 25216F) p 12 /16 ! 

~H 32768 + 174080P)^ 14 / 18 !-f... = °- ( 19 ) 

The minimum value of E seems to occur for about p 2 = 12, though it is 
not a sharply defined minimum ; however, with a similar approximation as 
in previous cases, we find the critical minimum of E to be 110. 

7. It is well known that, when fluid motion through a tube has changed 
from laminar to turbulent flow, the distribution of mean velocity over the 
cross-section alters so that the velocity becomes more nearly uniform over 
the greater part of the section, while falling to zero at the walls. This 
suggests a study of the comparative stability when the distribution of 
velocity alters in this manner, the boundary conditions being unchanged. 

However, it must be noted that we assume the distribution to be a steady 
state which has been acquired under a law of force, which may be deter- 
mined from the hydrodynamical equations, so as to give the required form 
for IT. 

A simple form, which illustrates the points in question, is 

U = (l + 1/2»)(1-^). (20) 

As n is made larger, the velocity approximates more closely to the mean 
velocity, tJ, over the greater part of the cross-section, while remaining zero at 
the walls. The corresponding law of force is, in the usual notation, 

X = v(4tn 2 -l)(TJ/a 2 ) v 2n - 2 . . (21) 

The greater the value of n, the more is the field of force concentrated near 
the walls, quite apart from the value of the viscosity. The flow approxi- 
mates to a uniform stream, but retaining the condition of zero velocity at 
the walls. 

The usual case of flow under a uniform field of force is given by n = 1. 

It is sufficient for comparison to work out another numerical case, say n = 2. 

We have then 

U = A(i-V). (22) 
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Equation (7) becomes 

(& -l]\jr~2]e(2cc^-i-3^) = 0, (23) 

\aor / \ dec / 

where ft = 5iH/8p 5 . 

The boundary conditions are 

^ = ; cl^jr/du = ; a — +j». 

Solving (23) by a power series 2iA n u n /n !, we have 

A* +6 = 2A n+4 -A n+ 2 + 2A(w4-l)(wH-2)(27H-3) A n . (24) 

As in the simpler cases, it is sufficient to choose fundamental solutions 
involving only even powers of a ; denoting these by i^ and i|r 2 we have 
f o = 1 + a 2 /2 ! + a 4 /4 ! + (l + 12ft) ofi/6 ! + (l + 192ft) a 8 /8 ! 

+ (l + 1032ft)a 10 /10! + (l+3552ft + 20160ft 2 ;a 12 /12 !•+•(! + 9492ft 

+ 696960ft 2 ) a 14 /14 ! + (1 + 21504ft + 8162256ft 2 ) « 16 /16!+... 
^ 2 = a 2 /2! + 2^/4! + 3aV6! + (4 + 168ft)a 8 /8! + (5 4-1656ft)a 10 /10! 
+ (6 -f 8184ft) a 12 /12!-h(7 + 28392ft-r 574560ft 2 ) a 14 /14! 
+ (8 + 78960ft + 11204352ft 2 ) a 16 /16 ! + (9 + 188496ft 

+ 102266496ft 2 )a 1 7l8 !-|-.... 
From the boundary condition 

Tfrodfa/doi — fad'tyo/doL = 0, 
we obtain the equation 

p + 2p*/3l + 8p 6 /5 !+32p 7 /7 ! + 128p 9 /9 ! + 512p u /ll ! 

+ (2048 + 129024ft 2 )j> 13 /13! + (8192 + 3280896ft 2 )p 15 /15 ! 

+ (32768 + 7753296ft 2 ) / 7 /17! + ... = 0. (25) 

Using this as an equation for E, we find by trial that the minimum value 
occurs near p 2 == 3 ; and the critical minimum value of E is 280 approximately. 

The corresponding value for the ordinary parabolic distribution (n = 1) is 
117. Thus, the critical value of E increases as the flow approximates more 
closely to a uniform stream, without slipping at the walls ; and, in this sense, 
the motion becomes increasingly stable. 

8. It has been stated that, under the boundary conditions u = 0, v = 0, 
there is thorough stability, in the ordinary sense, for simple shearing motion 
and probably also for laminar flow between fixed planes. In view of the 
behaviour of actual fluids in similar conditions, another suggestion has been, 
put forward by Hopf.* He proposes to express the influence of a wall by 
making the extra normal pressure, due to the disturbance, constant at the 

* L. Hopf, 'Ann. der Phys.,' vol. 59, p.- 538 (1919). 
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wall, together with no tangential slipping ; in fact, his boundary conditions 
come under case (ii) of § 4, namely u = 0, f yy = 0. With these assumptions, 
he applies the method of small vibrations to simple shearing motion between 
a fixed plane and a parallel moving plane. It appears that the motion is 
unstable for disturbances whose wave-length , exceeds a certain value; for 
smaller wave-lengths it is stable or unstable according to the value of E. 
Thus the motion is not thoroughly stable. Without discussing how far these 
assumptions express the behaviour of actual fluids and boundaries, we may 
see how they affect the energy method. 

We shall take the case of laminar flow between fixed planes, for which the 
previous calculations are available. 

The stream function -ty satisfies equation (11), and the boundary conditions 
are 

u = 0; —p + 2/jb?jvf(h/ = 0. 

From the equations (1), these are equivalent to 

n = 0; pvdU fdy—ficPu/dy 2 = 0, 

or, in the present notation, 

d^rjda = ; d*^ldaP--Zkaty = ; «= ±p. (26) 

Using the solutions ^ and yfr 2 , these give 

f o , (t2 y// ™2/^2)-^2 , (to ,,, -- 2/ ^fo) = 0, (27) 

where accents denote differentiation with respect to a. 

From the previous work, this equation involves odd powers of /.;. But & is 
3£R/4p 3 and w T e have to determine E in terms of p from (27). It follows that 
in this case there is no real solution of the problem of finding the critical 
minimum of E. 

It seems probable that it is only those motions which are completely stable 
in the ordinary theory which lead also to a real minimum for E. The suggestion 
may be stated in this manner : if a fluid motion is thoroughly stable when 
considered by the method of small vibrations applied to equation (2) and the 
boundary conditions, then it also possesses a real minimum value of E found 
from equation (7) and the boundary conditions. It has been pointed out 
that the latter equation is derived directly from the former, and it may be 
presumed that the minimum value of E depends in some manner upon the 
rates of decay of elementary vibrations and so may be used as a measure of 
the degree of stability of the system. 



